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Abstract. The one-dimensional Asymmetric Exclusion Process (ASEP) is a 
paradigm for nonequilibrium dynamics, in particular driven diffusive processes. It 
is usually considered in a canonical ensemble in which the number of sites is fixed. We 
observe that the (7ranc!-canonical partition function for the ASEP is remarkably simple. 
It allows a simple direct derivation of the asymptotics of the canonical normalization 
in various phases and of the correspondence with One- Transit Walks recently observed 
by Brak et al. 
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1. Introduction 

Driven pliysical systems tliat are out of equilibrium are characterized by steady states 
in which currents may flow. Examples of such systems are provided by the kinetics of 
biopolymerization p], transport across membranes j2] and traffic flow The existence 
of a macroscopic steady-state current has an important consequence for the microscopic 
modeling of such systems. Dynamics which are microscopically reversible (i.e. related 
to an energy function through the detailed balance condition) do not, by definition, 
admit macroscopic currents. Hence to understand nonequilibrium — and, in particular, 
driven — systems, it is useful to study closely microscopically irreversible processes. 

There is, at present, no unified statistical mechanical theory underlying 
nonequilibrium systems, and this presents a major barrier against progress towards 
a general, deep understanding thereof. Nevertheless, there is a class of one- dimensional 
driven-diffusive systems — typified by the Asymmetric Exclusion Process (ASEP) defined 
below — that have been exactly solved and shown to exhibit nontrivial steady-state 
phenomena such as phase transitions jl], spontaneous symmetry breaking [5^ and 
jamming [Hj. More recently, there have been some advances in extending notions from 
equilibrium statistical mechanics to nonequilibrium steady states. These include the 
Lee- Yang theory of partition-function zeros and phase transitions [3 IHl El smd, in a 
separate development, the derivation of nonequilibrium free energy functionals [TUj . 

It might be that the success of these approaches is due to the ASEP's steady state 
being related to certain types of random walks. This is because the latter are objects 
that are meaningful in the framework of equilibrium statistical physics — see, e.g., ^T] 
for a comprehensive overview. The connection between the steady state of the ASEP 
and One- Transit Walks was investigated in jT21 E| . In particular, it was shown that 
the transfer-matrix representation of the partition function can be interpreted in the 
framework of a similar formalism which proved instrumental for the exact solution of 
the ASEP 13]. In the present work, we show that the connection to random walks 
is transparent in an ensemble in which the system size of the ASEP fiuctuates. In the 
equilibrium picture, it is the lengths of the random walks that fiuctuate, and so we deem 
the ensemble grand canonical. In the course of this work we will show that the grand- 
canonical partition function has a surprisingly compact form, but one that nevertheless 
elegantly encodes the phase behavior of the ASEP in the canonical (fixed system size) 
ensemble. We conclude our work with some speculation as to the existence of more 
elementary solutions of the ASEP and related models based on the compactness of this 
formula. 

We begin by recalling the deceptively simple dynamics of the ASEP with open 
boundaries, indicated in Fig. 1. Particles are introduced at a rate a per unit time at 
the start of an A^-site chain and leave at a rate f3 at the other end. They can hop with 
probability one per unit time to the right if the space is empty, otherwise they remain 
stationary. Although the model is a nonequilibrium construct, defined purely in terms of 
microscopically irreversible transitions, it still displays steady-state solutions of various 
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Figure 1. Dynamics of the ASEP. The arrow labels indicate the rates at which the 
corresponding transitions occur; site labels are also indicated. 
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Figure 2. Phase diagram of the ASEP. 



sorts. A key quantity in describing the behavior of the model is the current of particles, 
which acts as an order parameter in the different phases. 

The model displays a phase diagram that possesses both first- and second-order 
phase transitions. In Fig. 2, phase (i) is a high-density phase with a density profile that 
decays exponentially towards the right boundary and has a current J = (3{1 — f3). Phase 
(ii) on the other hand is a low-density phase that has an exponential decay in the density 
profile from the left boundary and a current J = a{l — a). The line a = /5 < | is a first- 
order transition line at which the current exhibits a discontinuity in its first derivative. 
Finally in phase (iii), the current assumes a maximal constant value of J = |. Since J 
has a discontinuity in its second derivative passing from phase (iii) to either (i) or (ii) 
the horizontal and vertical boundaries are second-order transition lines. 

It is at present unknown whether one can define for general nonequilibrium steady 
states a partition function from which all thermodynamic information can be obtained. 
Nevertheless it is known that a quantity constructed by summing the statistical weights 
(unnormalized probabilities) f{C) over all configurations C in the steady state of the 
ASEP does at least capture the model's phase behavior (by applying, e.g., a Lee- Yang 
analysis jH]). Denoting this analog of the partition function by Z, we have 

Z = J2f{C), (1) 

c 

where the weights f{C) are defined by the microscopic transition rates W{C ^ C) in 
the model and the requirement that probability fluxes in the steady state balance: 

[f{C')W{C' ^ C) - fiC)W{C ^ C')] = . (2) 

The probability of a configuration C is then given by P{C) = f{C)/Z. Since the weights, 
f{C), are fixed only up to an overall scale by (j2I), one can ensure that Z is uniquely 
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defined by, for example, insisting tliat it is polynomial in the transition rates W{C C) 
with any factors common to all weights f{C) removed I13j. 

The tour de force in |lj gives exact formulae for these weights for an ASEP with N 
sites. The resulting expression for the canonical partition function, Z^v, reads 

^ m{N-p)\ (l//3)-(l/a) • ^ ' 

The current is given by J = Ztv-i/^at and various correlation functions can also be 
obtained [3]. For later convenience we note that combinatorial factor appearing in Q 
is given by the Ballot numbers B^^p 

_ p{2N-l-p)\ 
^^'^ - N\{N-p)\ ' 

so that 

^"-^"""^^ (1//?) -(!/«) 

N P 1 1 

-yBN.y-^- (5) 

p=l q=0 

Ballot numbers arise in the enumeration of Dyck paths which are walks on a 
diagonally-rotated square lattice with the following properties: (a) they are constrained 
to lie above or on the x-axis; (b) they can only move in the north-east or south-east 
directions; and (c) they start and end on the x-axis. The quantity Bj^^p is the number 
of Dyck paths with length 2N and p returns to the x-axis, including the final return 
that is always present pLlj . 

With this in mind, (0) can be interpreted as a weighted sum over paths as follows. 
Take a Dyck path from (0, 0) to (2A^, 0) with p returns, and from the q^^ return 
onwards (g = 0,1,..., p) reflect the path about the x-axis. The resulting construct 
is a One-Transit Walk, an example of which is shown in Fig. El The weight, a~'^ (3~^p~'^\ 
associated with the walk can be gleaned by assigning a factor 1/zi = 1/a to each 
contact with the x-axis from above, excluding an initial upward step, and similarly a 
factor 1/^2 = 1//? to each contact from below, excluding a final upward step. Note that 
initial or final downward steps do contribute factors of 1/^2 and 1/2:1 respectively since 
these correspond to the cases q = and q = p- 

A slightly different scheme is that implemented in |12j. There, a fugacity, zi, is 
assigned to each down step and Z2 to each up step, apart from those ending on the 
X-axis. The partition function of this model is 

Af p 

which has the same form as Zn up to the analytic factor of {ziZ2)^ ■ This prefactor, 
which brings Q into the standard form proposed in ^3 E], has no bearing on the 
critical behavior. 
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Figure 3. A One- Transit Walk obtained by composing two Dyck paths. Each contact 
with the X-axis from above, apart from the start, gives a factor 1/a and each contact 
from below, apart from the end, gives a factor of 



2. A Grand Canonical ASEP 

We now switch to an alternative viewpoint of the ASEP which allows both the phase 
behavior and the connection to One- Transit Walks to be seen in a very straightforward 
way. The key step is to relax the constraint of a fixed system size N. To this end we 
introduce a fugacity, z, which is conjugate to and define a grand-canonical partition 
function, Z, as 

oo 

Z{z) = J2 Znz"" . (7) 

After substituting in the expression (jH)) for Zn and changing the order of summation, 
one can use 15^ 

2z ) p\{N + p)\ ' 

to give the following remarkably simple result for Z: 

[x{z) — a){x{z) — p) 
In this expression x{z) is the generating function for Catalan numbers. 



xiz) = 2 (1 - VT^) , (10) 

which crops up in numerous combinatorial contexts ^H] , such as counting well-bracketed 
words or bifurcating trees. Writing Z{z) explicitly in terms of z, we have 

Z(,) = ^ (U) 

(l-2a-v/T^4^)(l-2/5-v/r^4i) ' ^ ' 

or, when a = P, 

One can alternatively use the Lagrange inversion formula with x = z/{l — x) to show 
that expanding Z{x{z)) also gives the desired result in ((Tj). 
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We remark that the expression (jllj) has recently been derived in two distinct 
contexts. First, it turns out to be instrumental in determining the joint distribution 
of density and current in the ASEP fTl]. To obtain the grand-canonical expression, 
a novel "relaxed" operator algebra method was invoked |18|. Meanwhile, in a version 
of the ASEP in which the dynamics allow the system size to fluctuate (thereby giving 
physical meaning to the grand-canonical ensemble), an expression similar to, but more 
general than |TT|) arises by determining a quantity related to the resolvent of a matrix 

Standard results for the asymptotics of expansions of functions [20] allow us to 
read off the known large N behavior of Zn |1] directly from (fTT|) and ()12|1 . The large 
behavior of is dominated by whichever pole, a; ~ a or x ~ /5, is closest to the 
origin, or, if a, /3 > 1/2, by the square root singularity in x{z) itself. The first two cases 
correspond to regions (i) and (ii) in the phase diagram (Fig. 2) and the latter to region 
(iii). In the latter case, with a < P for example, the square root singularity in a/1 — Az 
in (fTTj) gives rise to a 4^/A^^/^ asymptotic behavior and substituting in z = 1/4 picks 
up the coefficient 

1 1 



4^ a(3 



_(2«-l)2 (2/5-l)2j ' ^^^^ 
in perfect agreement with |1] . The a = (3 result in region (iii) follows in a similar fashion 
from 

In regions (i) and (ii), on the other hand, it is the pole terms which lead to the 
dominant asymptotic behavior. For (ii), the pole at x = jS, z = i3{l — j3) is dominant 
since a > (3 and the square root singularity at z = 1/4 has not yet been reached. In 
this l/(/5(l — (3))^ behavior results from the pole term, giving 

Finally, on the second-order boundary lines (a = 1/2, P > 1/2) and (/? = 1/2, a > 1/2) 
the singularity is of the form 1 / \/T^^^ and gives an asymptotic behavior of the form 
4Af/]vi/2_ ^ 1/2, (3 > 1/2), for instance, 

4^ 2/3 

~ vrV2iVi/2 2p - 1 • (^5) 
There is a close link between the grand-canonical partition function for the ASEP 
and the random walk picture, as we now show. The generating function for Dyck walks 

with a fugacity w for each contact with the x-axis is known to be given by [TT] 

w 

g"'"'^)- 1-^(1- VI -4.)/2 - ^'''^ 

in which ^1/2 is the fugacity for each step. Comparison with ()1H) reveals that Z can be 
written as 

Z{z) = ai3 GB(l/a, z) Gb(1//3, z) . (17) 

The product form of this expression implies that the grand-canonical ensemble of paths 
comprises two weakly-interacting Dyck paths of fluctuating length, one with a contact 



The Grand- Canonical Asymmetric Exclusion Process and the One-Transit Walk 



7 



fugacity w = 1/a and the other with contact fugacity w = One way in which 

such paths can be reahsed is as One- Transit Walks as shown in Fig. El Then, the a/3 
prefactor simply takes account of the unweighted initial and final steps. Furthermore, 
from ()17|). it is clear that the critical behavior of the ASEP is a combination of 
that for two adsorbing Dyck walks where the second-order lines correspond to the 
adsorbtion/desorbtion transition for one of the walks and the first-order line to a co- 
operative transition not seen for a single walk. These results are in perfect accordance 
with those recently derived in for the canonical ensemble. 



3. Comparison with parallel dynamics 

It is interesting to compare the grand-canonical generating function in (jllll with that 
describing an ASEP with parallel dynamics, which was derived in using recursion 
relations rather than a direct summation 



aP{l + pz) 2(1 - p){ap - p^z) - al3h'^{l - pz) - a^h'^ ^ {I + pzf - Az 
ZJz) = — 



2p^{l - [3){1 - a){z - zm){z - Zld) 



;i8) 



where 



and 



^[(l-p)-(l-a)(l-/5)] (19) 



a(p — a) , , 

Zld = -if^ ^ 20 

P'^{1 — a) 



_ Pip - P) 

- p^il-P) ^^^> 
are the poles which govern the low density and high density phases respectively. 
With a parallel update all possible particle moves into empty spaces are carried out 
simultaneously with probability p at each time step, rather than the random sequential 
updates implicity used elsewhere in the discussion. The parallel update ASEP is a 
special case of the Nagel-Schreckenberg [221 niodel for traffic flow and the transitions in 
the model can be interpreted as jamming transitions for the traffic in a one dimensional 
flow. The phase diagram has a similar topology, and similar phases, to the random 
sequential update ASEP |2I1I2S15 but the second order lines now lie at a, /3 = 1 — ^^1 — p. 

It is clear from the definition of the dynamics that substituting a = pa, P = p P 
and taking the limit p — > will recover random sequential sequential updates. If we 
take this limit in (jlSp we find 



Z„^o{z) = ap- 



l-a-p + 2ap + {!-&- py\-Az 



+ 0{p) (22) 

{z - a{\ - a)){z - p{\ - P)) ^ ' ' 

where the leading term is just a rationalized form of the grand-canonical generating 
function derived earlier in (jlip . The limiting process for obtaining random sequential 
dynamics from parallel dynamics can thus be carried out at the level of the grand- 
canonical generating function. 
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It is also noteworthy that it is possible to write 




(23) 



where X-{z) is the negative root of 

x{p — x) 
— x) 



(24) 



The x{z) appearing in the grand-canonical generating function for random sequential 
dynamics is one of the roots of z = x{l — x) , so the structure of the parallel grand- 
canonical generating function above is similar. In fact, it turns out that the lattice paths 
described by this latter generating function generalise the Dyck walks of the random- 
sequential case in a natural way [21] . 

4. Discussion 

Taking stock, we have found that the introduction of an initially ad-hoc grand-canonical 
approach to the ASEP leads to an extremely simple, unified description of the entirety of 
the model's steady-state phase behavior, not to mention a straightforward confirmation 
of lattice path equivalences derived for the canonical ensemble in |12| IT!^ . It is, in 
fact, possible to go further and relate properties of the lattice paths — obtained using 
standard equilibrium statistical mechanical techniques — to observables in the ASEP's 
nonequilibrium steady state. For example, density profiles JH] and correlation functions 
have been calculated in this way. Hence it is evident that the lattice path 
interpretation of a nonequilibrium partition function can be developed into a useful 
calculational tool. 

In this work we employed the grand-canonical ensemble as a mathematical device. It 
is, in fact, possible to give physical meaning to this ensemble if one introduces processes 
that allow the size of the system to fluctuate at a rate controlled by the fugacity z. A 
model with such processes was investigated in In that work, the A^-site chain is 
thought of as a 'cluster' embedded in a larger system of particles that can attach to and 
detach from the cluster at its boundaries at various rates. For a particular choice of these 
rates, in which z is the ratio of the attachment to detachment rate, the normalisation 
given in corresponds to our expression (|TT|l . If particles attach at a rate greater than 
the current of particles through the cluster, the size of the cluster explodes to infinity, 
a phenomenon that explains the divergence of the grand-canonical partition function at 
a critical value of z equal to the internal current. In ^Hj, an analogy was made to a 
gas-liquid transition. 

Given the compactness of the grand-canonical normalisations for the ASEP with 
both random-sequential (fTTj) and parallel dynamics (j2Sl), one might suggest that 
the grand-canonical approach might be of general utility in establishing equivalences 
between nonequilibrium and equilibrium steady states. One way in which this hypothesis 
might be given further weight would be to derive the analogous quantity for the Partially 
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Asymmetric Exclusion Process (PASEP) I77| |2H] in which hops to the left at a rate 
q are also entertained. This model contains an additional phase transition at g = 1 
which is of a different nature to the boundary-driven transitions shown in the ASEP 
phase diagram, Fig. 2. It may well be that a grand-canonical partition function for the 
PASEP would shed additional light onto this transition. 

Finally we remark that the simple form of (fTT|) . which required previously known 
results as input, is somewhat tantalizing. In particular it hints at the possibility that it 
might be directly calculable from purely physical principles. This would be an avenue 
worth pursuing since, as noted in the introduction, partition functions appear to have 
some relevance to nonequilibrium steady states, yet there are as yet no systematic 
approaches for obtaining them. 
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